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We demonstrate how the matrix-valued mode-coupling theory of the glass transition and glassy 
dynamics in planar confinement converges to the corresponding theory for two-dimensional (2D) 
planar and the three-dimensional bulk liquid, provided the wall potential satisfies certain condi¬ 
tions. Since the mode-coupling theory relies on the static properties as input, the emergence of 
a homogeneous limit for the matrix-valued intermediate scattering functions is directly connected 
to the convergence of the corresponding static quantities to their conventional counterparts. We 
show that the 2D limit is more subtle than the bulk limit, in particular, the in-planar dynamics 
decouples from the motion perpendicular to the walls. We investigate the frozen-in parts of the 
intermediate scattering function in the glass state and find that the limits time t —¥ oo and effective 
wall separation L — ► 0 do not commute due to the mutual coupling of the residual transversal and 
lateral force kernels. 
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I. INTRODUCTION 

Confined liquids have been studied extensively in 
physics, in particular their phase behavior [1, 2], dynam¬ 
ical properties [3- 13], and their structural characteriza¬ 
tion [14-16]. Furthermore confined liquids are interme¬ 
diate between a low-dimensional and bulk liquid and dis¬ 
play an intriguing interplay of near-range local ordering 
and the confining length. In case of a slit geometry, i.e., 
two parallel flat hard walls with effective separation L 
one can elaborate the full range from two-dimensional 
(2D) up to three-dimensional (3D) liquid behavior by 
varying the wall separation from zero to infinity. For 
instance, for equilibrium phase transitions one can study 
the crossover [17] from the Kosterlitz-Thouless transi¬ 
tion [18-23] in 2D liquids to the conventional phase tran¬ 
sitions of a 3D liquid. Similarly, the crossover behavior of 
the glassy dynamics of a quasi-two-dimensional confined 
liquid towards a bulk system provides insight in the na¬ 
ture of the mechanism of structural arrest. 

In bulk systems many features of the slowing down of 
structural relaxation upon cooling or compression have 
been rationalized in terms of the mode-coupling theory 
of the glass transition (MCT) developed by Gotze and 
co-workers [24, 25]. The predictions of MCT include the 
emergence of a non-trivial long-time limit of the interme¬ 
diate scattering functions, called glass-form factors asso¬ 
ciated with two-time fractals in the vicinity of the criti¬ 
cal point. Particularly, MCT entails a sharp dynamical 
glass transition ,e.g., for a hard-sphere liquid at a criti¬ 
cal packing fraction. Although in nature this transition 
is smeared, the various predictions of MCT in three di¬ 
mensions have been confirmed in the past two decades 
by experiments and computer simulations [25, 26]. 

In order to study the dependence of the glass transition 
on the spatial dimension MCT has also been worked out 
for two-dimensional single-component [27] and binary liq¬ 


uids [28-30]. Qualitative [27, 29, 30] and quantitative [30] 
comparisons of the MCT results for these 2D liquids with 
those from experiments [31] and simulations [30] have 
been reported. The two-dimensional systems are found 
to behave qualitatively similar to their three-dimensional 
counterparts. 

Motivated by the numerous experimental and simu- 
lational results for the glass transition of confined liq¬ 
uids [3, 6-13], mostly for a slit geometry, MCT has been 
generalized recently to describe dense liquids in planar 
confinement [32-34]. In contrast to the MCT for the 
glassy dynamics in bulk or disordered structures [35, 36], 
the inhomogeneous packing in the slit requires to consider 
symmetry-adapted matrix-valued intermediate scatter¬ 
ing functions to characterize the density fluctuations of 
the confined liquid. In the case of hard spheres, the tran¬ 
sition line, separating the regime of collective frozen-in 
states from liquid states, has been determined as a func¬ 
tion of the slit width. An intriguing multiple reentrant 
transition for wall separations on the scale of a few parti¬ 
cle diameters has been predicted and corroborated by 
recent molecular-dynamics simulations with controlled 
polydispersity [37]. 

While MCT has been successfully tested for planar, 
bulk and confined liquids, the natural question arises: 
Does MCT for liquids in planar confinement for L —> 0 
and L —> oo properly converge to MCT for 2D and 3D liq¬ 
uids, respectively? Particularly, one would like to know 
if the dynamical behavior of a strongly confined liquid 
is approximately described by a two-dimensional system. 
To provide an answer is the major goal of the present 
work. In particular, we will demonstrate how the matrix¬ 
valued MCT due to the inhomogeneous structure of the 
confined liquid reduces to the MCT of scalar intermediate 
scattering functions for the two-dimensional [27] and the 
bulk case [25]. In contrast to the bulk limit, the planar 
limit L —> 0 turns out to be rather subtle. 
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The outline of this work is as follows. In the next sec¬ 
tion we introduce the model, the quantities of interest, 
and recall the equations of motion for confined liquids 
within the MCT approximation. In Sec. Ill we discuss the 
behavior of the MCT functional for small wall separation 
in terms of the proper convergence of the static proper¬ 
ties towards their two-dimensional counterparts demon¬ 
strated recently [38]. In particular, we show that the 
planar MCT is recovered for all finite times in the limit 
of vanishing plate separation. Next we study the fixed- 
point equation for the glass-form factors for small slit 
widths. In Sec. IV we demonstrate that the MCT equa¬ 
tions include the bulk behavior as limiting case as the 
wall separation becomes infinitely large. Section V pro¬ 
vides a critical assessment of the different convergences 
and possible implications for the glassy dynamics in ex¬ 
treme confinement. 


II. CONFINED LIQUIDS: BASIC QUANTITIES 
AND MCT 


The microscopic setup, the derivations of the equations 
of motion of the relevant dynamic quantities of interest is 
detailed in Ref. [33]; here we summarize its main features, 
in order to keep the present paper self-contained. 

Consider a liquid of N particles of mass m between two 
parallel, planar walls with cross section A and separation 

H. Then a point in phase space is specified by the set of 
coordinates parallel and perpendicular to the wall x n = 
(f n , z n ) and corresponding momenta p n = ( P n , P*), n = 

I, ..., N. We use impenetrable walls with additional wall 
potential of the form U({z n }]L) = Yln=i M( z n', L) and 


U(z-,L) 


fW w (z) for|^| < L/ 2, 
|oo for| 2 :| > L/2. 


(1) 


Here, the effective wall separation L is introduced as the 
transverse length accessible to the particles, and therefore 
we distinguish between point particles and hard spheres, 


L = 


H — a, 
H, 


hard spheres, 
point particles. 


( 2 ) 


The interaction energy of a particle with either of the 
walls U w {z) is assumed to be smooth, with possible di¬ 
vergencies at z = ±L/2. In principle, both walls can 
interact differently with the liquid, i.e., asymmetric wall 
potentials U w {z) = U-(L/2 + z) + U + {L/2 — z) with 
IA _ (a:) yf U+ ( x ) are allowed. Periodic boundary condi¬ 
tions are imposed parallel to the walls in x-y direction. 
We suppress the parametric dependence on the wall sepa¬ 
ration L in the following to allow for a compact notation. 
The pair interactions V(x) = V(r, z) depend only on the 
mutual distance |x|, i.e., V({x n }) = J/n<m V(|x„ -x m \). 
Then the Hamilton function is specified by 


H({p n },{xn}) = Y 


N r -a 
Pn 


n=l 


.2 TO 


+ 1A(z n ) 


Throughout this paper the microscopic dynamics are as¬ 
sumed to be Newtonian generated by the Hamiltonian 
H({Pn}, {x n }), hi particular, collisions with the flat hard 
walls are elastic thereby conserving momentum parallel 
to the walls. The trajectory in the N particle phase 
space is denoted by {{r n {t)},{z n (t)},{P n (t)},{P^(t)}) 
and all calculations are performed in the canonical en¬ 
semble. The thermodynamic limit N —>• oo, A —> oo is 
anticipated for fixed 2D number density no = N/A. With 
the accessible volume of the particles V = AL , we find 
for the 3D number density n = N/V = no/L. 

The modulation of the equilibrium density profile n(z) 
in the slit is encoded in discrete Fourier components 

= J dzexp(iQ /i z)n( 2 ), (4) 

where the mode index /j refers to discrete wave numbers 
= 27r/r/L, p £ Z and integration is performed over 
the accessible slit 2 G [—L/2, L/2]. 

The fundamental variable of interest is the microscopic 
fluctuating density mode 

N 

Pu.(q,t ) = ^exp[iQ^„(f)]e 19 ' r " (t) , (5) 

n =1 

where q = ( q x , q y ) are the conventional discrete (for finite 
A) wave vectors in the x-y- plane. Particle number con¬ 
servation relates the time derivative of t ) to the cur¬ 
rent densities parallel and perpendicular to the walls [33] : 

1 N „ _ 

(£ t ) = — Y &Q (q ■ p n (t), Pn (t )) exp [i QfiZn (t )] e' 9 ' 1 '" {t) , 

TYl 

n— 1 

( 6 ) 

with channel index a =||, J _Here, the short-hand no¬ 

tation for the unit vector q = q/q and the selector 
b a (x,z) = x5 a \\ + zS a _i has been employed. The emer¬ 
gence of several decay channels is reminiscent of the 
mode-coupling theory of molecules [39, 40], where gen¬ 
eralized density modes couple to both translational and 
rotational currents. 

The basic quantity of the MCT of liquids in slit geom¬ 
etry is the generalization of the intermediate scattering 
function 

Sn»(q,t) = 0 )). ( 7 ) 

We shall make use of a natural matrix notation 
[S(q,t)]^ = ( q , t ), and similarly for other correla¬ 

tion functions throughout this paper. MCT requires 
the Fourier coefficients n M , the static structure factors 
S l j, ll (q) = S liv {q,t = 0) and the static current-current 
correlators 


[J(q)C = ^{q) 


= 

= 6 afl v 


n 0 


+ V({x n }). (3) 


( 8 ) 
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as known input [33] with thermal velocity vth = 
(/csT/m) 1 / 2 . General properties of the static and dy¬ 
namic correlators including symmetry relations have 
been discussed in detail [33]. 

Employing the Zwanzig-Mori projection operator for¬ 
malism [25, 41] exact equations of motion for the collec¬ 
tive correlators S llw (q. t) can be derived [32, 33]. For 
later purposes it is more convenient to consider the 
Laplace-transformed equations, where the convention 
for the transformed matrix correlators is S^ v (q,z) = 
if™ 0 dte lzt S^(q,t), Imfs] > 0 1 . The equations in 
the Laplace domain for S^ v (q,z) deal with general¬ 
ized matrix-valued fraction representations. First, we 
express S fJjl ,(q,z) in terms of current-memory kernels 
K[iv(q,z) [33], 

S(q,z) = - [^S _1 (g) + S^(g)K(g,2)S -:l (g)] , (9) 

which split by the perpendicular and parallel current to 

k^{q lZ )= Y. b a (q,Q,)IC^(q,z)k(q,Q„). (10) 

a/3=||,-L 

The kernel [ JC(q , z)]™& = ( q , z) can be represented by 

K{q, z) = ~ [zJ^^q) + J~ 1 {q)^XH{q, z)J~ 1 {q)] , 

( 11 ) 

which involves the force kernel z). The con¬ 

stitutive MCT ansatz expresses the force kernel in the 
time domain in terms of the intermediate scattering func¬ 
tions [33], 

wf„(q,t) = ^3 

X (q 1 , t)S ^ V2 (g 2 , t)X ^ UlV2 (q, qiq 2 )*- 

( 12 ) 

The vertices have been calculated relying on a systematic 
convolution approximation [33], 

= — Nv^-^S^q!+q 2 [6 (<7 ■ <?1, QfJ,— M 2 ) 

xc ^-<i 2 ,w(9i) + (1^2)], (13) 

and involve the direct correlation functions c ^ (q) , which 
are related to the static structure factor S lw (q) by the 
inhomogeneous Ornstein-Zernike equation [34, 42]. 

The Eqs. (9)-(13) involve the initial conditions 
Snv(q,t = 0) = S^(q) and K.ffl{q,t = 0) = J$(q) 
thereby constituting a complete set of coupled non-linear 


1 The frequency z should not be confused with the transversal 

variable z. The distinction is clear from the context. 


equations which have to be solved self-consistently. The 
investigations performed in the following rely on the 
proofs that the mode-coupling equations provide the ex¬ 
istence of unique solutions, which has been demonstrated 
rigorously for Brownian dynamics in monocomponent 
simple liquids [25, 43] and mixtures [44] and only re¬ 
cently for the MCT with multiple relaxation channels, 
where the here discussed confined MCT is merely a spe¬ 
cial case [45]. 

III. CONVERGENCE TO THE PLANAR MCT 

In this section we discuss how the MCT equations con¬ 
verge towards the MCT for 2D liquids, which is one of 
the major results of this work. Let us emphasize that 
this requires to discuss the theory for small but finite 
slit width L > 0. Convergence implies that fluids con¬ 
fined to small slit widths behave similarly to the truly 
two-dimensional case, the errors can be made arbitrarily 
small upon decreasing the slit width. 

In the limiting regime of small wall separation the in¬ 
plane motion is presumed to be close to the dynamics of 
a two-dimensional system. Furthermore, of all the struc¬ 
tural properties entering the MCT equations the two- 
dimensional structure factor should dominate the equa¬ 
tions for small slit widths. The MCT equations encode 
the confinement by the walls via the structural input in 
terms of the density modes and the symmetry-adapted 
static structure factors both as initial values for the cur¬ 
rent correlators and for the intermediate scattering func¬ 
tions and on the associated direct correlation functions 
via the vertices in the MCT functional. One subtlety 
arises since for small wall separation L > 0 the wave 
numbers Q M = 27r/r/L associated with the perpendicular 
motion blow up. 


A. Density profile and static correlators 

Since the MCT relies on the structural properties of 
the fluid as known input, the question of convergence to 
a planar limit is intimately related to the structure of 
fluids in extreme confinement and the rapidity of the ap¬ 
proach to a two-dimensional system. This issue has been 
addressed in two of our earlier works [38, 46] in terms 
of a cluster and cumulant expansion in case of hard core 
and smooth potentials, respectively. There, one can not 
only estimate the order of convergence but actually calcu¬ 
late the leading corrections to both thermodynamic and 
structural quantities with respect to a planar reference 
system. Here we shall be interested only in the leading 
terms and quantify the order of the corrections in terms 
of Landau symbols O(-) and o(-) as the slit width be¬ 
comes small L —> 0. 

The convergence of the structural quantities has been 
demonstrated [38] assuming analytic wall potentials ful¬ 
filling the smoothness criterion U{z = zL)—U( 0) = O(L), 
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which states that the particle-wall interaction should be 
controlled for fixed scaled transverse coordinates z as the 
slit width approaches zero. Then the density profile be¬ 
comes flat [38] even on the scale of the plate distance 
n(z) = [no/L)[l + 0(L)]. This in turn implies conver¬ 
gence for the Fourier modes of the density, 


n p = 


no = const. 
O(L) 


for n = 0, 
else. 


(14) 


Note, that the Fourier mode no is independent of L from 
the sum rule n 0 = f d zn(z) = N/A. We assume here 
the worst case of asymmetric walls, whereas for symmet¬ 
ric walls 0(L) can be replaced by the faster convergence 
0{L 2 ), see Ref. [38]. Let us repeat a word of warning 
here. The existence of the convergence is not guaran¬ 
teed in general and, in particular, potentials diverging at 
the walls, e.g., Lennard-Jones and Coulomb potentials, 
do not belong to the class indicated above. Yet, pure 
hard walls trivially fulfill the smoothness criterion and 
they constitute the reference case we wish to address. 
From the convergence properties of the density modes, 
Eq. (14), one can immediately infer the convergence of 
the static current correlator, Eq. (8), to 


Jffin) = S' 



for /i = iz, 
else, 


(15) 


which becomes to leading order diagonal with respect to 
the discrete modes. 

Similarly, for the static structure factors it has been 
shown [38] that 


S^{q) 


5(g)[1 + 0(L 2 )} for fi = v = 0, 
(1 — + @(L) else, 


where S(q) denotes the static structure factor of the cor¬ 
responding 2D liquid. Thus the structure factor matrix 
[S(g)] M „ = S l _ lv (q) becomes diagonal as well in the 2D 
limit and in addition, S /llJ ,(q) for /r ^ 0 becomes ideal- 
gas-like for L —» 0. 

The direct correlation function of the confined liquid 
in the slit converges to 

Cn V {q) = L 2 [c{q)5 ll . 0 S v0 + rfc^iq) + o(L 2 )], (17) 

where c{q) is the corresponding 2D direct correlation 
function and the correction amplitude c M „(g) is inde¬ 
pendent of the slit width L. The prime observation of 
Ref. [38] was that the corrections are 0(L 2 ) irrespective 
of the wall potential, whereas for the structure factors for 
(n,v) 7^ (0,0) or the density profile the leading correc¬ 
tions are O(L). This latter property plays an important 
role for the convergence of MCT for confined liquids. 


B. Mode-coupling theory: t finite, L —> 0 

We start the investigation of the convergence of the 
confined MCT towards the planar MCT by discussing 


the memory kernels. We assume first, that times t 
and frequencies z are held fixed independent of the wall 
separation, while the limit L —» 0 is performed. For 
this purpose, it is favorable to introduce rescaled modes 
Q M = LQ M = O(L 0 ), which are independent of the slit 
size. Using the convergence of the direct correlation func¬ 
tion, Eq. (17), neglecting terms of 0(L 2 ) from Eq. (13), 
one obtains 

/u/i2 (?’ *?1 ! 92) 

= - Nv th n ofiq,qi+q 2 {fi a Hq ' 9l)c(gi)6 Ml 0'W 

+ LS a± Q l _ l - l _ l2 c fJj - fJj2tl _ ll (qi) + (lo 2) j. (18) 

Thereby, we find that the leading contribution of the ver¬ 
tices in the limit L —> 0 stems from a =||. Keeping now 
only the leading order, i.e., we set L = 0, and assuming 
a priori that S l _ lv (q,t) = O(L 0 ) for all {n,v), in the 2D 
thermodynamic limit N —> 00, A —> 00 with no = N/A 
fixed, the nonvanishing memory kernel assumes the fol¬ 
lowing form 

{q,t) =y v th J j^{[q-qic(qi)] 2 Soo{qiA)S^(q 2l t) 

+ [q ■ qic{qi)][q ■ q 2 c(q 2 )]So v {qu t) 

+ (lfi 2)}, (19) 

where q 2 = q — qi- In particular, for /i = v = 0 the 
memory kernel contains only couplings of Soo (q, t) 

/ j 2 

qiq 2 )Soo(qi,t)Soo(q 2 ,t), (20) 

and the vertices V(q, qiq 2 ) coincide with the ones of the 
two-dimensional theory [27] 

V(q,qi,q 2 ) = Wi)c(gi) + (1 O 2)] 2 . (21) 

The remaining memory kernels assuming S^ ( q , t) = 
O(L 0 ) for all (fi, v) are of higher order in L, viz. 


mi 1 (9, t) 

mm^) 

= 0(L), 
= 0(L ), 
= 0{L 2 ) 

as one infers from Eq. (18). 



The notable property of the MCT functional, 
Eqs. (19), is that if the intermediate scattering functions 
Snv(q,t) are diagonal in the mode indices, this property 

is preserved by SJljil (q, t). Since also the static structure 
factors and current correlators, which serve as initial con¬ 
ditions, are diagonal to lowest order in L the equations 
of motion, Eqs. (9) —(11) do not generate off-diagonal 
terms. More formally one can show that all time deriva¬ 
tives d z 5 Ml/ (g, t)/d< z | t=0 , l € N 0 are diagonal, similar to 
Ref. [47]. Since the solutions have been demonstrated 
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to be unique [45], the thus constructed solution remains 
diagonal for all times t > 0. 

In particular, one finds that the equations of motion 
for Soo (q, t) decouple completely from the remaining di¬ 
agonal ones as we demonstrate below. To simplify no¬ 
tation we drop the mode indices and write S(q,t) = 
Soo(q,t) and similarly for the static structure factor 
S(q) = <Sbo(<Z, t = 0). Furthermore we make contact with 
the notation of Ref. [27]: 


^ooO?^) = toqVt h m(q,t), (23) 


with the characteristic frequency Q 2 = q 2 v 2 h /S(q). Us¬ 
ing this result and Eq. (15) it follows from the second 
fraction representation, Eq. (11), that its solutions are 
diagonal in (a,/3) and in (/x, i/), and that the equations 
for K,^l(q,z) and K.^(q,z) decouple for all /x. Conse¬ 
quently, one obtains a closed equation for £jj]j(<7, 2) in¬ 
volving rh(q,z ), only. Abbreviating the 2D relaxation 

kernel K(q,z) = Koo(q,z) = q 2 K,H(q,z) one finds from 
Eq. (11) 


K(q,z) 


2 2 
9 «th 


z + Q 2 rh(q , z) 


(24) 


Taking advantage of the diagonality of S IJ/V (q) and Eq. (9) 
one obtains the well-known double-fraction representa¬ 
tion 


S(q,z) 


_ -S{q) _ 

~ - 02/ [ z + 02 m(q,z)\ ' 


(25) 


Going back to the temporal domain leads to the gener¬ 
alized harmonic oscillator equation 


S(q, t)+ttqS(q, t)+Og f m(q,t—t')S(q,t')dt' = 0, (26) 
Jo 


coinciding with the MCT equation of a 2D liquid [27]. 

It is interesting to ask what the MCT equations yield 
for the remaining diagonal correlators S^^^q, t) for /x ^ 0 
for small plate separation L. While the equations for 
Soo(q, t) allow for a direct limit L = 0, the L-dependence 
cannot totally be eliminated for the other quantities. We 
evaluate Eq. (11) again for L = 0, and the contraction 
with the selectors, Eq. (10), implies for the current ker¬ 
nel K^^q, z) = —Qf l L~ 2 v 2 h /z. The transversal coherent 
scattering function then becomes 


Snn(q,z) = 


- 1 


z-Q\L- 2 v 2 Jz 


for /x ^ 0, (27) 


which is the equation of motion of uncoupled undamped 
oscillations. Thus for L —>• 0 and /i / 0 the correla¬ 
tors S „„ ( q , t) display fast harmonic oscillations with fre¬ 
quency = Q tll L~ 1 v th := 47r/x/r(T,T) and are inde¬ 
pendent of the planar wave number q. Their basic pe¬ 
riod r = t(T,L) = 2L/vt\i is just the time a single par¬ 
ticle with thermal velocity t> t h needs for a single bounce 


between both walls. This time scale becomes small for 
L —> 0, thus it is separated from the microscopic dynam¬ 
ics r <C to = 1 /fig of the planar dynamics. Indeed, if one 
inspects Eq. (25) on the time scale r, i.e., zt = 0(1), the 
planar correlator S(q,t) for t/r = 0(1) has not evolved 
for L —> 0: S(q,t) = S(q). On this rapid scale the MCT 
kernels are negligible, which is in accordance with the no¬ 
tion that they are designed to describe the slow dynami¬ 
cal processes leading to structural arrest. If one performs 
the same reasoning as described above neglecting the 
memory kernels for the MCT equations of motion of the 
tagged-particle correlator in Ref. [34], then one obtains 
the same equation of motion for the transversal dynam¬ 
ics S${q, z) = -1/ [z - Q] l L~ 2 v 2 Jz\ for /x ± 0, if 
the tagged particle is of the same species as the host- 
liquid particles. Thus, the fast transversal dynamics of 
the collective intermediate scattering function reduces to 
the dynamics of the incoherent scattering function, i.e., 
Sjj*il(q,z) = (g, 2:) for /x ^ 0 and L —> 0 indicating 

that the residual perpendicular dynamics is merely a one- 
particle dynamics. 

We conclude that the equations of motion are capa¬ 
ble to account for a dynamical decoupling of lateral and 
perpendicular degrees of freedom, splitting off the pla¬ 
nar glassy dynamics from a fast one-dimensional ideal- 
gas like motion in a finite box. Nevertheless, we expect 
that the true transversal dynamics is not correctly con¬ 
tained in the MCT equations, but the fundamental time 
scales are still properly reflected. 


C. Mode-coupling theory: L finite, t — ¥ oo 

The present case allows to study the glass form factors 
as a function of the wall separation L. Beyond a cer¬ 
tain critical point, the MCT equations describe structural 
arrest characterized by nonvanishing glass-form factors, 
i.e., for confined liquids F^(q) := lim^oo S^q, t) ± 0. 
In experiments and simulations, these frozen-in parts de¬ 
scribe the plateau values of the intermediate scattering 
function in the dense or supercooled regime [26, 48, 49]. 
We inspect the solutions for the glass-form factors in the 
limit L —>■ 0. From Eq. (26) one can readily extract 
the limit lim^oo lim^o Soo{q, t) = F(q), which coin¬ 
cides with the glass-form factors of the planar MCT [27] 
in case of structural arrest. Interchanging the limits, 
i.e., taking first the limit t —> oo and then L —» 0, 
liiui-j-o lim t _).oo >S/u'(fo t) = lim L ^. 0 F fJjV (q), these glass- 
form factors differ qualitatively from those obtained from 
MCT in 2D [27], as will be demonstrated in the following. 

The fixed-point equation for glass-form factors F^ v (q) 
have been derived in Ref. [33] by performing the limit 
z —>• 0. Since frozen-in force kernels display a pole at 
zero frequency, Tl^(q,z) = ~{^{q)/z)[l + o{ 1)] for 
z —> 0 one obtains 

F (q) = [S” 1 ^ + S- 1 (g)N- 1 (g)S- 1 ( 9 )] _1 , 


(28) 
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with contractions 

[ivV(5) = E &“(<?» £ -1 6 i.) 

«/3=||,J- 

x[*r- 1 {q)]'*bP{q,L-'Q v ), (29) 

and the inverse of the frozen-in part of the MCT func¬ 
tional 

[A r-\q)]f v := (30) 

It has been proven that these equations exhibit a unique 
maximal solution [33, 45]. 

The goal here is to find a self-consistent solution for 
the glass-form factors F /lu (q) in terms of estimates in 
powers of L. The strategy is to perform a conver¬ 
gent iteration suggested in Ref. [33] and to keep track 
of the respective orders in L. We initialize the iter¬ 
ation with F /liy (q) = O(L 0 ) for all /z, v. Then, with 
the vertices, Eq. (18), one infers for the long-time lim¬ 
its limt-Kxj (q, t) = $^{q), that the mode-coupling 

functionals display orders 3j[|(g) = O(L 0 ), (?) = 

0(L), ^(g) = 0(L) and ^j~(q) = 0(L 2 ). In con¬ 
trast to the frequency-dependent equations of motion, 
Eq. (11), the limit L = 0 cannot be performed, since 
then the force-kernel matrix $(q) becomes singular and 
inversion in Eq. (30) is not possible. 

Keeping the leading-order estimates for the kernels one 
obtains immediately by Eqs. (28) —(30) the estimates 
F 00 (q) = O(L 0 ), F 0v (q) = 0(L 2 ) { V ± 0), FM = 
0(L 2 ) (/i ^ 0) and F^(q) = 0(L 4 ) 0) as first 

iterate. Reinserting the first iterate in the mode-coupling 
functional does not reduce the orders further, as demon¬ 
strated in detail in Appendix. Hence the solutions, 


Foo(q) = O(L 0 ), 



(31) 

F,o(q) = 0(L 2 ) 

for 

M 7^ 0, 


F 0v (q) = C(L 2 ) 

for 

v + 0, 


Ffiv(q) = 0{L 4 ) 

for 

M,*' 7^ 0, 



are the unique solutions characterizing the glass states 
in extreme confinement for small but finite L. In fact 
the Landau symbols O(-) can be replaced by asymptotic 
proportionality ~ (•), i.e., the orders cannot be improved. 

Let us emphasize, that for all estimates lim^^o F^(q), 
the frozen-in parts of the force kernels (q) for 
a, [3 = || _L are mutually coupled. Thus, the equations for 
t —> oo first do not decouple for small wall separations, 
in striking contrast to the case where the limit L —>• 0 is 
performed for fixed finite times. In particular, this ob¬ 
servation entails for the glass-form factor hin^^o Foo{q), 
that it does not coincide with the glass-form factors of 
the two-dimensional MCT, Eq. (26); see the discussion 
in Appendix below Eq. (A8). 

In conclusion, we have demonstrated, that the equa¬ 
tions of motion display a delicate dependence in the lim¬ 
its L —> 0 and t —>• oo. Both limits do not commute, 


i.e., lim^oo lim^^o (q, t) ± lim L ^ 0 lhn t _>oo S^{q, t). 
This result necessarily implies the existence of an In¬ 
dependent diverging time scale tl, on which the lateral 
dynamics couples to the transversal one. We emphasize 
that the iteration scheme for the nonergodicity parame¬ 
ter, see Ref. [33, 45], remains valid for arbitrarily small 
L > 0 and yields always the solutions, where lateral and 
transversal degrees of freedom are coupled. 

IV. CONVERGENCE TO THE 
THREE-DIMENSIONAL BULK 

In this subsection, we demonstrate that the MCT for 
confined fluid approaches the standard MCT of the 3D 
glass transition for wall separations approach infinity pro¬ 
vided the static correlations are short-ranged. We discuss 
first the convergence of the static structure factors and 
the density profile. Then, we use these properties to ex¬ 
tract the bulk limit of the MCT for confined liquids. 

A. Density profile and correlators 

We assume that the fluid-fluid and fluid-wall interac¬ 
tions are short-ranged and show that the bulk behavior 
for the average density and the two-particle static correla¬ 
tion functions is attained as the wall separation becomes 
large. This limit has to be performed such that the 3D 
density n — no/L remains fixed. The spatial dependence 
of the density n(z) arises due to the wall potential U{z) 
and the pair-potential V(\x — x'\). It displays a significant 
variation only in the vicinity of the wall, and becomes 
constant otherwise. Therefore we can write 

n(z) = n + An(z ), (32) 

where A n(z) decays on a scale of a wall correlation length 
£ away from the wall, which can be defined by 

£ = — [ dz|An(z)|. (33) 

n J 

We recall that n is the bulk density. In general £ depends 
on L , but for large wall separation L —> oo it assumes a 
finite limit depending on temperature and density, only. 
We assume that the confined liquid is sufficiently far away 
from a critical point. Then we obtain with Eq. (4) and 
no = nL 

nii =n 0 [6 M o + 0{^/L )\, (34) 

since the integral [cf. Eq. (4)] over A n(z) is of order £. 

Similar considerations apply for the distinct part 
G( d ) z , z') of the density-density correlation function 
G(r,z,z’) (see Ref. [38] for conventions). We write the 
distinct part as asymptotic expansion with respect to the 
corresponding bulk correlator 

G (d) (r, z, z') = i [G3 D (z) + A G^ (r, z, *')] , (35) 
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where x = ( f,z ) and G^(x) denotes the bulk dis¬ 
tinct part of the density-density correlation function [50] , 
which displays rotational and translational invariance. 
The corrections AG^(r, z, z') again decay on a wall cor¬ 
relation length which we assume to be of same order as 

e 

The Fourier decomposition 2 of the distinct and the 
corresponding self part [G^ (r, z, z') = n(z)S(r)S(z — 
z')/n 0 \ follows by Eqs. (34) and (35) 

S$(q) = SW(k)6^ + 0(Z/L), 

S${q) = S^ + O^/L). (36) 

Here we adopt the convention that k = (q. Q ^) abbrevi¬ 
ates the 3D wave vector, whereas q is reserved for two- 
dimensional vectors, in particular k = \k\ = (q 2 + Q 2 ) 1 / 2 . 
For the structure factor we arrive at 

SM = s(k) [S„ v + 0(£/L )}, (37) 

where S(k) is the structure factor of the bulk liquid. 
From the Ornstein-Zernike relation in confined geome¬ 
try [34, 42] one immediately obtains 

c lu/ (q) = L[c(k)6 lu/ + 0^/L)], (38) 

with the bulk direct correlation function nc{k ) = 1 — 
1/S(k) [50]. Thus, the static structure factors S^{q) and 
the direct correlation matrix (c Mi ,(g)) become diagonal 
for L —> oo. 


B. Mode-coupling theory 

For the bulk limit of the MCT equations of motion we 
also require the convergence of the static current corre¬ 
lator (q). From Eqs. (8) and (34) we obtain 

(9) = + 0(Z/L ). (39) 

Thus all static correlators become diagonal with respect 
to the discrete mode indices /r and v in the bulk limit L —> 
oo. Substituting c^(q) from Eq. (38) into the memory 
kernel [Eq. (12)] and using Eq. (39) we find 

W°P(q,t) =^ n Vh 

qi Ml 1 '! 


X 

b a 

7 Q/Xi 

) + (1 O 2) 



V q y 


X 

» 3 

K ‘ ,Q^i)c{qi,Q^ 

+ (1 2) 


X Vl (l^l, t^S^ 2L , 2 ((^2 ; ^)j (40) 


2 The convention for correlators in slit geometry is A llv (q) = 
Jd 2 rdzdz r A(r, z, z') exp(— iQ^z) exp {iQ v z')e~ l x r 


where q 2 = q — qi,/i 2 = /t — n\, V 2 = v — v\. The mem¬ 
ory kernel appears to be nondiagonal both in (a, /3) and 
(/r, v). However, as in the two-dimensional limit we shall 
show that the subspace in which S lll/ (q,t) is diagonal, 
i.e., 


Sp V (q,t) = S(k,t)5p V (41) 


remains invariant under the MCT equations. 

The assumption, Eq. (41), is consistent with the diago- 
nality of the static structure factor S^ v (q) [Eq. (37)], the 
initial condition for S/j, v (q,t). One infers for the memory 
kernel the simplification 


hE 


k i 

g- 9l 
9 

9'9i 
9 


; Q/J. 1 J c (k l) + (1 H 2 ) 
, Qfi i) c(fci) + (1 -H- 2 ) 


x S(ki,t)S(k 2 ,t)Sf_, 


(42) 


with = k — k\ and the corresponding wave vectors read 

ki — (gz, Q 

The thermodynamic limit L — > oo ,A -A oo,N —> oo 
such that N/AL = n restores isotropy in addition to 
homogeneity. Therefore we can choose the direction of 
the ’external’ wave vector k = ( q ., Q M ) in a suitable way. 
Without restricting generality we choose k = (q = 0, Q^), 
which implies for the projections q ■ q% = 0 . With k = 
\k\ = \Q^\ and Q w = k ■ ki, the selector simplifies to 


b a (q-qi,Q fJ . i ) = k-k i 5 a - L . (43) 


Then the only nonvanishing elements of the memory ker¬ 
nel in the thermodynamic limit become 


9^(9, *) = <W 


/ 


with the 3D vertices, 


d 3 fci 


V(k, kik 2 )S(ki,t)S(k2, t), 

(44) 


V{k,kik 2 ) = -u t 4 h (k ■ fci)c(fci) + (1 2) 


-i 2 


(45) 


Again we make contact to established notation and write 
DJl^(q,t) = S^nlv^m^kA) with flj? = k 2 v 2 h /S(k) and 
m(k,t ) coincides with the MCT kernel for a 3D liq¬ 
uid [25]. 

The remaining steps are similar to the 2D case. Re¬ 
placing in Eq. (43) Q by = fc, the total current 
correlators reduces to K^(q = 0 ,z) = k 2 K,jq^(q = 0 ,z) 
[see Eq. (10)]. By Eqs. (39) and (11) it becomes diagonal 
Kqviq = 0, z) = M f 2 h <5 M ,, in the discrete wave numbers. 
We abbreviate K(k,z ) = K^iq = 0,z) = k 2 /C^(q = 
0, z) and arrive at the equation of motion 


K(q,z) 


"th 


z + f2|m(fc, z) 


(46) 



















Diagonality of S^(q) and of K^(q, z ) in Eq. (9) implies 
that S /liy (q = 0, z) = S(k,z)5 /ll , remains diagonal, i.e., 
the assumption of Eq. (41) remains consistent. One ob¬ 
tains again a double-fraction representation 


S(k,z) 


~S{k) 

z - n 2 k / [z + n 2 k rh(k, z)Y 


(47) 


Finally, the Laplace back transfrom of Eq. (47) yields the 
3D generalized harmonic oscillator equation 


S(k,t) + Cl1S(k,t) + f m(k, t — t')S(k, t')dt' = 0, 

J o 

(48) 

which is the well-known MCT equation for glassy dynam¬ 
ics of a bulk liquid [25]. Since for large wall separations 
we expect S^ v (q,t) = 0(1;/L) for /j ^ we will not 
discuss the MCT equations for these correlators which 
describe the dynamics close to the walls. In contrast to 
the 2D limit, the limits t —> oo and L —> oo do commute, 
since the fixed-point equation for the nonergodicity pa¬ 
rameter, Eqs. (28)—(30), converges properly for L —> oo 
to its bulk counterpart [25]. 


V. SUMMARY AND CONCLUSIONS 

We have demonstrated that the mode-coupling equa¬ 
tions for confined liquids [32, 33] for finite times con¬ 
verge to the MCT of the planar and bulk liquid for ef¬ 
fective wall separation L —> 0 and L —> oo, respectively. 
In particular, we have shown how these limiting cases 
emerge from the reduction of the matrix-valued theory, 
which accounts for the inhomogeneous structure within 
the slit. In both cases, time is initially finite, while the 
limits L —> 0 and L —> oo are performed, respectively. 
The recovery of the correct limits demonstrates the con¬ 
sistency and robustness of the MCT ansatz also for con¬ 
fined liquids. 

Several conditions are imposed on the static level to 
guarantee the existence of these limits. For the 3D 
limit we have assumed short-ranged particle-particle and 
particle-wall interactions and the thermodynamic state 
of the liquid is located sufficiently far away from a crit¬ 
ical point. For the 2D limit, we have required that the 
particle-wall interaction is such that structural properties 
approach their two-dimensional counterparts, in partic¬ 
ular, the density profile becomes flat sufficiently fast for 
L —> 0. The analyticity of the particle-wall interaction is 
sufficient for this condition. 

The generalized matrix-valued intermediate scattering 
functions become diagonal in both limits. For finite times 
in the limit L —> 0 we have demonstrated that the 2D 
glassy dynamics of the lateral degrees of freedom decou¬ 
ples for L —>• 0 from the transversal dynamics. While 
the dynamics parallel to the walls coincides with the 2D 
MCT dynamics, the latter display fast dynamics on a 
time scale L/v th mimicking ideal-gas-like motion in a one¬ 
dimensional box of size L. On this short time scale we 


find that the mode-coupling contributions to the trans¬ 
verse force fluctuations can be neglected for the descrip¬ 
tion of this rapid dynamics. It becomes obvious that the 
transversal collective scattering function in this regime 
reduces to the incoherent scattering function indicating 
that this dynamics is a mere one-particle dynamics. More 
generally, we expect that the MCT equations for the in¬ 
coherent dynamics (see Ref. [34]) also converges to its 
planar and bulk case displaying similar features as found 
for the collective correlators discussed here. 

The structural properties and thermodynamic phase 
behavior of extremely confined fluids can be determined 
using an effective two-dimensional pair potential, as has 
been shown recently [38, 46]. Hence it would be inter¬ 
esting to use the two-dimensional MCT with these ef¬ 
fective potentials and to compare to the corresponding 
results of the MCT in confinement. Similarly, one could 
compare computer simulations in slit geometry with two- 
dimensional simulations using the effective pair poten¬ 
tial. An approach in the same spirit has been pursued 
for colloid-polymer mixtures where the effective Asakura- 
Oosawa depletion interaction has been employed also 
for studying dynamical properties. Similar strategies 
have been applied also for asymmetric colloidal mix¬ 
tures [51, 52] or star polymers [53]. We anticipate that 
this works reasonably well also in the present context 
although a microscopic justification is lacking. 

Let us discuss the limit of small wall separations in 
more detail, since the planar limit displays peculiarities. 
Interestingly, we find that the limits t —> oo and L —^ 0 
do not commute. If the limit L —> 0 is performed first 
for fixed time, then applying the limit t —> oo leads to 
the glass-form factor of the corresponding 2D glass state. 
Taking the reverse order yields a glass state, where resid¬ 
ual mutual couplings between the perpendicular and par¬ 
allel frozen-in stresses remain. This subtle dependence 
on the performed limits suggests the existence of a diver¬ 
gent L-dependent time scale tl, on which the transver¬ 
sal ideal-gas-like motion couples to the lateral degrees of 
freedom. This time scale occurs as a consequence of the 
non-commutativity of the limits t —> oo and L —> 0 of 
the underlying MCT equations in confined geometry and 
provides an interesting prediction for future simulations 
and experiments. 

Let us speculate on the relevance of such a time scale 
tl diverging for L —> 0 which competes with the struc¬ 
tural a-relaxation time. For instance, it becomes conceiv¬ 
able that for sufficiently small L and fine-tuning of the 
density, the lateral degrees of freedom are frozen up to 
times t tl and behave effectively as a two-dimensional 
glass. At later times the coupling to the transversal de¬ 
grees of freedom sets in and yields either a different glass 
state or even melts the glass entirely. Hence the transver¬ 
sal fluctuations effectively soften the planar interaction. 
Such a (partial) melting scenario differs from glass-glass 
transitions or reentrant transitions, e.g., for attractive 
glasses [52, 54, 55] or binary mixtures [56, 57], since it 
is a purely dynamical phenomenon for a single thermo- 
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dynamic state. The emergence of two different regimes 
of the glassy dynamics is based on the assumption that 
the transversal degrees of freedom are slow variables even 
in extreme confinement. Therefore, it would be interest¬ 
ing to test this prediction experimentally or by computer 
simulation and verify that two different regimes, sepa¬ 
rated by a divergent time scale tl, indeed behave quali¬ 
tatively differently. Potential candidates for experimen¬ 
tal realizations are superparamagnetic colloidal particles 
at a liquid surface [31, 58-60] where small fluctuations 


in out-of-plane directions yield a generic mechanism to 
weakly couple to transversal degrees of freedom. 
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Appendix A: Glass-form factors in extreme confinement 

In this Appendix we investigate the convergence of the glass-form factors, i.e., lirii£_>o F^ v (q). We demonstrate here 
that the estimate F 00 (q) = O(L 0 ), F 0v {q) = 0(L 2 ) (: v ^ 0), F^q) = 0(L 2 ) (/j, ^ 0) and F^(q) = 0(L 4 ) (n,v ± 0) 
yields a consistent solution for the fixed-point equation; see Eqs. (28) —(30). The required ingredient for this analysis 
is the asymptotic behavior of the vertices for L -A 0, Eq. (18). For all a,/3, we list here the expanded functionals to 
lowest required powers in L [convention q 2 = q — <fi]: 
a = p=\\: 


stlo?) [(q-qi) 2 c{qi) 2 F 00 (qi)F^{q 2 ) 

Qi 

+ (q- 5i )(q- q 2 )c(qi)c(q 2 )F 0v (q 1 )F fJ . 0 (q 2 ) + (1^2) 

+ [(^■qi) 2 c(qi)K-u 2 ,u 1 (qi) F o^(qi)Fu t u 2 (q2) 

qi vi v-2 

+ (q ’ <h)(q ' q 2 )c(qi)c*_ Vl V2 (q 2 )F 0l , 1 (q 1 )F lll , 2 (q 2 ) + (1 -(A 2) 

4 

+ l2 2N^2 [(9-9i) 2c (?i)V/*2,/ii(9i) f B o(«i)^2492) 

qi M 1 M 2 

+ {q- qi)(q- Q 2 )c{qi) ,M2 {q2)Fp 1 o{qi)F li2l/ (q 2 ) + (1 tA 2)J, 

« =\\,P =-L: 

$lt{q) qi) c (qi)Q^-^K- V2 , Vl (qi) F o V i(qi)F^ 2 {q2) 

qi V1IS2 

+ (q2) F ou 1 {qi ) F ^ 2 (<? 2 ) 

qi V1V2 

+ (1 -o- 2), 

a =_L, /3 =||: 


(Al) 


(A2) 


qi M1M2 
v A 

(q- 92)c(Q , 2)0 A 6- ( u 2 c M _ Al2lAll {qi)F lilV {q{)F tl2 Q{q 2 ) 

qi MlA»2 

+ (1 1A 2). 


(A3) 
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a = /3 =_L: 


v 4 

S^j, (?) = L 2 jy ^ / Qn—itaQv—v-J'li— M 2 ,Mi (9l)^i/—i/a,i/i (9 i)-Pmi"i (9l)-f/i 2 ‘ / 2 fe) 


<?i 


+ i‘ 


7 ; 4 
2 _th 


2N EE Q fx — [l 2 Qh' — l'l &H — M2,Ml (^l)^l/—1/1,1/2 (92)^11/1 (gi)i^2,2 (92) 
+ ( 1 « 2 ). 


Qi Ml M2 
VI M2 


(A4) 


We insert our estimates for the glass-form factors, Eq. (31), and count the leading orders in L. It turns out that 
for estimating the orders, it is sufficient to distinguish between mode index zero (’O’) and a generic non-zero mode 
(’O’). In this short-hand matrix notation, we find 


W(«)) = 



II 0 

II 0 

± 0 

T 0 

II 0 

' O(L 0 ) 

0(L 2 ) 

0(L 3 ) 

0(L 3 )1 

II 0 

0(L 2 ) 

o\l 4 ) 

0(L 5 ) 

o{l 3 ) 

T 0 

0(L 3 ) 

0(L 5 ) 

0(L 6 ) 

o\l 4 ) 

T 0 

L 0{L 4 ) 

0(L 3 ) 

0(L 4 ) 

o\l 2 ) _ 


(AS) 


The inversion of this matrix can be done by introducing block matrices 


W(9)) 


_L 


A 

C 


± 

B 

D \ ’ 


e-g-, 


etc. Block-matrix inversion yields 


([ r-\q)]^) = 


A- = (A^y) 


0 0 
0 ro(L°) o(L 2 y 
0 \_o\l 2 ) o[l 4 ) J’ 


II 

_L 


II 

(A-BD-'C)- 1 

D 4 C(A - BD 1 C)~ 1 


± 

A~ 1 B(D CA^B)- 1 
(.D - CA^B )- 1 


(A6) 


(A7) 


(AS) 


An important insight is, that the estimates of each entry, e.g., compare A and BD~ 1 C , are of the same order in 
L. The consequence is, that each entry of ((S’” 1 (9)]“^) contains contributions from all force kernels Suf (?) with the 
various combinations of a and (3. In particular, we emphasize that for the entry (a = /3 =||; /z = v = 0) the scalar 

(BD 4 C) 00 = £^[(S^)- 1 ]^' = (A9) 

K-y 


is of the same order as Aoo = Sod, the latter corresponding to the arrested 2D MCT kernel. This in turn implies that 
Aqo — (BD~ 1 C )00 does not coincide with the expected 2D result: 


Aoo - {BD~ l CU y $jjj}. 


(A10) 


aa-W) (All) 



II0 

II 0 

± 0 

± 0 

II 0 

' O(L 0 ) 

0(L~ 2 ) 

0(L~ 3 ) 

0(L-y 

II 5 

0(L~ 2 ) 

0(L ~ 4 ) 

0(L~ 5 ) 

0(L~ 3 ) 

T 0 

0(L~ 3 ) 

0(L~ 5 ) 

0(L~ 6 ) 

0(L” 4 ) 

T 0 

ML- 1 ) 

0(L~ 3 ) 

0(L ~ 4 ) 

0(L~ 2 )- 


Explicit inversion leads to 
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From the diagonality property of the current correlator, Eq. (15), ([A/" (?)]u5) retains the same estimates [see 

Eq. (30)]. Contraction with the selectors, Eq. (29), yields in matrix notation, 


U iV o [0(L~ 2 ) 0(L~ 4 )_' 


(A12) 


Since the static structure factor is diagonal to leading order, see Eq. (16), one infers from Eq. (28) the final result: 


0 0 
0[ O(L 0 ) 0(1?) 
- Q | o(L 2 ) 0(L 4 ) 


(A13) 


where we have recovered our initial ansatz. Since the solutions are unique [45], we have found the consistent solution 
for L —> 0. In particular, F 00 (q) couples to all residual kernels of S'^(g) [see explicitly Eqs. (A9) and (A10)] and 
thereby the such obtained glass-form factor differs from a purely two-dimensional glass state [27] . 
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